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ABSTRACT 

We measure the photometric properties of 105 giant arcs that were identified in systematic searches 
for galaxy-cluster-scale strong lenses in the Second Red-Sequence Cluster Survey (RCS-2) and the 
Sloan Digital Sky Survey (SDSS). The cluster lenses span 0.2 < z; < 1.2 in redshift, with a median 
zi = 0.58. Using broadband color criteria we sort the entire arc sample into redshift bins based on 
u — g and g — r colors, and also r — z colors for the ^ 90% of arcs that have z— band data. This analysis 
yields broad redshift constraints with 71^4 % of the arcs at z > 1.0, 64I4 % at z > 1.4, 56^4 % at 

z > 1.9, and 21^2 % at z > 2.7. The remaining 291°^ % have z < 1. The inferred median redshift 
is Zs = 2.0 ± 0.1, in good agreement with a previous determination from a smaller sample of brighter 
arcs {g < 22.5). This agreement confirms that Zs = 2.0 ± 0.1 is the typical redshift for giant arcs 
with 5 < 24 that are produced by cluster-scale strong lenses, and that there is no evidence for strong 
evolution in the redshift distribution of arcs over a wide range of band magnitudes (20 < g < 24). 
Establishing that half of all giant arcs are at z > 2 contributes significantly toward relieving the tension 
between the number of arcs observed and the number expected in a ACDM cosmology, but there is 
considerable evidence to suggest that a discrepancy persists. Additionally, this work confirms that 
forthcoming large samples of giant arcs will supply the observational community with many magnified 
galaxies at z > 2. 

Subject headings: gravitational lensing — galaxies: high-redshift galaxies — cosmology: observations 
— large-scale structure of universe 



1. INTRODUCTION 

The statistics of strong gravitational lensing by 
clusters of galaxies provides a test of cosmological 
models, by comparing the observed abundance of giant 
arcs against predictions from ray-tracing cosmologi- 
cal N-body simulations. Incidents of strong lensing 
are often identifiable by the formation of giant arcs, 
which are background sources that are strongly lensed 
into multiple, often merging, images by a foreground 
gravitational potential. The frequency of giant arcs is 
a complex observable that scales with the halo mass 
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Based on observations taken at the Southern Astrophysi- 
cal Research Telescope (SOAR), a collaboration between CNP- 
Brazil, NOAO, The University of North Carolina at Chapel Hill, 
and Michigan State University and the Canada-France-Hawaii 
Telescope (CFHT) on Mauna Kea, which is operated by the Na- 
tional Research Council of Canada, the Institut National des Sci- 
ences de rUnivers of the Centre National de la Recherche Scien- 
tifique of France, and the University of Hawaii. Additional sup- 
porting observations come from the 2.5m Nordic Optical Tele- 
scope, operated on the island of La Palma jointly by Denmark, 
Finland, Iceland, Norway, and Sweden, in the Spanish Observa- 
torio del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias, and the Gemini Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: The United States, The United Kingdom, 
Canada, Chile, Australia, Brazil and Argentina 



function, the detailed properties of high-redshift galaxy 
populatio n, and the interna l properties of the halos 
( e.g., iGrossman fc NaravanI [l98 8: Bartcl mann et all 
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Bartelmann et al.l ()1998[ ) first suggested that the num- 
ber of giants arcs on the sky is under-predicted by 
ACDM by approximately an order of magnitude, and 
later observational comparisons using small samples of 
giant arcs corroborated the apparent over- abundance 
of gia nt arcs on the sky compared to ACDM predic- 
tions (iLuppino et alj [19991: iZaritsky fc Gonzaleg 120031: 
iGladders et all 120031 : iLi et all 120061 1 In principle, it 
should be possible to produce predictions for giant arc 
statistics that match what we observe, assuming that 
we possess a suitably accurate model describing the 
formation and evolution of structure over cosmic time. 
In practice, the apparent failure of predictions for giant 
arc statistics to match the observations by as much as 
an order of magnitude is clear evidence that more work 
remains to be done. 

Resolving the persistent discrepancy between the 
abundance of giant arcs observed and the number pre- 
dicted by early simulation efforts requires progress on 
two fronts: 1) improved empirical constraints from much 
larger samples of giant arcs that are identified in uniform 
data with a well-characterized selection, and 2) refined 
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predictions based on simulated giant arc samples. It is 
tempting to consider the possibility of giant arc statistics 
providing direct observational evidence fo r divergence 
from the concordance ACDM co smology (iFedeli et alJ 
l2008HD'Aloisio fc Natar aian"2011). but we must first ex- 
plore other physical explanations for the discrepancy be- 
tween giant arc c ounts in observations and simulations. 
Motivated by the iBartelmann et al.l ()1998[ ) result, theo- 
retical work over the past decade has explored a variety 
of effects which, if unaccounted for in theoretical pre- 
dictions, can result in an under-abundance of simulated 
arcs compared to the number of giant arcs observed in 
the real universe. 

Suggested effects include the contribution of baryons 
to cluster lensing cross-sections in the f orm of cen- 
tral massive galaxies and subst r ucture (iFlores et al.1 
200^, 'Mcncghctti ct al.'l2000|, [200l IHennawi et al.ll2007l: 
Mencghctti ct_aL 201Qi) . the effects of dark matter 
being dragge d into the cores of clusters by cool- 
ing baryons (iPuchwein et al.l 120051: iRozo et al.l 120081: 
iWambsganss et ah 20081) . triaxiality of cluster mass pro- 
files CQguri et al."2003': 'Dalai ct al.il2004l: IHennawi et al.l 
Il007; Mcncghctti ct al. 2010), accounting for addi- 
tional but uncorrelated structure in the universe 
that is projected alon g the hne-o f- sight of lens-source 
systems dW ambsganss et all l2005t iHilbert et all 120071: 
iPuchweiiifc Hilbertji2009l) . short time-scale increases in 
the strong lensing cross-section due to ma.jor merg- 
ers jlTorri et al.l 120041: iFedeh et al.l 120061: IHennawi etahl 
120071 ) ■ and the properties of the assumed background 
galaxy source populatio n ()Hamana fc Futamasd 119971: 
IWambsganss et ah! I2004D . The aforementioned papers 
have generally quantified the impact of various effects 
in the context of comparing predictions and obser- 
vations for the frequency of the formation of giant 
arcs, with several authors m aking arguments to dis - 
clai m the discrepancy noted by IBartelmann et alJ (1 19981) 
(e.g.. IWambsganss et al.ll200llHoresh et al.ll2005[|201l[ ). 
These arguments, however, are based on the identifica- 
tion of possible mechanisms for explaining or alleviating 
the apparent giant arc discrepancy, and there is a per- 
sistent lack of empirical evidence to confirm or deny the 
relevance of these different mechanisms. 

1.1. The Background Source Population 

In this paper we focus on the background galaxy pop- 
ulation as an input into predictions for giant arc abun- 
dances in simulations, and make a direct measurements 
of the giant arc redshift distribution for a large sam- 
ple of real giant arcs. E vidence from theoretical work 
(jWambsganss et all I2004D suggests that the use of ap- 
propriate distribution of background galaxies may have 
the greatest potential to explain the dramatic discrep- 
ancy between the giant arc counts observed and those 
predicted by simulations in a ACDM cosmology. The 
redshift distribution of background sources impacts the 
global efficiency for giant arc production through the an- 
gular diameter distan ce term in the gravitational lens 
equation (Equ. 13. ,Naravan fc Bartelma nn 1996), which 
has an explicit dependence on the ratio of the angular 
diameter distances between the lens and the source, and 
the observer and the source. In practice, this means that 
the critical surface mass density for a given foreground 
lens is a function of the source redshift, where higher 



source redshifts result in lower critical surface mass den- 
sities. When analyzing the total cross section for strong 
lensing of all clusters within some cosmological volume, 
a background source population that is shifted to higher 
redshifts will therefore require lower surface mass densi- 
ties of foreground structures in order to become super- 
critical (i.e. to be strong lenses). There has, however, 
been little effort made to systematically measure the 
background source redshift distrib ution of a well - defined 
sample of arcs until very recently. iBavliss et al.l (|2011a[ l 
first measured the redshift distribution for a spectroscop- 
ically complete sample of 28 giant arcs that were observed 
with Gemini-South, and they demonstrate that the giant 
arc redshift distribution can only be reliably measured 
using a large s ample of unif ormly selected giant arcs. 

[Wambsgans s et al.l (|2004[ ) note an increase in the to- 
tal strong lensing cross-section for galaxy cluster-scale 
lenses of a factor of ~ 10 — 20 when the background 
sources that are lensed into giant arcs are assumed to 
be at corresponding source planes of Zg = 1.5 — 2 rather 
than Zg = 1. Other work in which this effect is quantified 
finds smaller increases in the strong lensing cross-section 
by a factor of ~ 3 or so (Dalai et al. 2004; Li ct al. 2003 
iFedeli et al.ll2006( ). but this still implies a significant po- 
tential boost in the total cross-section for giant arc pro- 
duction if giant arcs are typically formed by lensing of 
galaxies at z. ^ = 2, rather than t he z,, = 1 source popula- 
tion used by IBartelmann et al.l (fl998). The importance 
of understanding the properties of the background source 
population that is lensed into arcs is further highlighted 
by more recent simulation work to generate predictions 
for giant arc statistics. 

In agreement w i th early semi- a nalyt ic arguments 
(lOg uri et al']l2003[ ). IHennawi et"all 12007) identify the 
dominant uncertainty in their ability to robustly predict 
giant arc counts to be the poor constraints on the number 
density of galaxies at high redshift that are available to 
be strongly lensed by foreground cluster lenses. This un- 
certainty is fundamentally a refiection of our ignorance of 
the faint tail of the surface brightness function of galaxies 
at high redshifts. When we co nsider the imphcation s of 
IWambsganss et al.. (200l ) and IHennawi et al.l (|2007D to- 
gether it becomes clear that we cannot realistically make 
any useful statements about cosmology from tests of gi- 
ant arc statistics without a firm understanding of the 
population of galaxies that are typically lensed into gi- 
ant arcs. 

The first spectroscopic measurement of the redshift dis- 
tribution for a complete sample of well-selec ted giant arcs 
indica tes a median arc redshift of Zg = 1.82 (|Bavliss et al.l 
l2011al) . However, the sample used in that work was still 
relatively small (Na^cs = 28), and potentially biased to- 
ward high-redshift giant arcs due to a target selection 
criteria which favored larger giant arc radii (Rare) sys- 
tems, where Rare is the average angular distance between 
an arc and the centroid of the foreground lensing poten- 
tial. There should be a net correlation between Rare and 
the intrinsic Einstein radius, R^;, of a given lens, and 
for a given lens the Einstein radius grows monotonically 
with the redshift of the lensed source. However, there 
should also be a large scatter in the correlation between 
Rare and He that could easily be washed out by other 
factors such as the shape of the critical curves for each in- 
dividual lens. It is not clear to what degree this potential 
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selectio n bias should impact the results of iBavliss et alj 
(|2011af) . but there is a somewhat reasonable argument to 
be made that the re dshift distribution measurement in 
IBavliss et alJ (|2011al ) could be biased significantly high. 
In order to ensure a robust understanding of the prop- 
erties of the galaxies that are typically lensed into gi- 
ant arcs, it is important that we extend these investiga- 
tions further to incorporate very large sample of arcs that 
are not subject to th e same selection effects as those in 
IBavliss et alJ (|2011a[ ). However, measuring spectroscopic 
redshifts for a sample of ~ 100 giant arcs is extremely 
observationally expensive, especially in a scenario where 
spectroscopic completeness (i.e. redshift measurements 
for ^ 100% of the observed sample) is important. 

We also note that counting giant arcs may be in 
large part free from the issue of magnification bias that 
is known to affect, for example, counts of unresolved 
galaxies or quasars behind fore ground lens i ng stru ctures 
(jBroadhurst et all I1995D . As lOalal et"all (|2Q04[ ) point 
out, the magnification due to gravitational lensing is spa- 
tial in nature, with surface brightness a conserved quan- 
tity. The magnification does not, therefore, impact the 
number density of spatially resolved objects on the sky. 
Giant arcs are always, by definition, spatially resolved 
along at least one (and sometimes both) axis(es) on the 
sky. Assuming that magnification bias does not strongly 
influence counts of giant arcs, we make no attempt to 
model or correct for the affect in this paper. In support 
of this assumption we note that § 14.2! includes a direct 
comparison of the median redshifts of two samples of 
giant arcs that were selected from imaging data that dif- 
ferent in depth by 1.5 — 2 magnitudes. The observed 
similarity in the median redshifts of these two samples 
provides some empirical evidence to support the argu- 
ment that magnification bias does not strongly impact 
surveys of giant arcs, but future arc searches that probe 
to fainter magnitudes will be necessary to robustly quan- 
tify the precise effect of the magnification bias on giant 
arc counts. 

1.2. Measuring the Properties of A New Giant Arc 
Sample 

The "giant arc problem" posed by the apparent excess 
of giant arcs on the sky relative to ACDM predictions 
persists because of limitations on both the theoretical 
predictions and observational constraints. Making se- 
rious progress toward addressing the issue requires im- 
proving the fidelity of the machinery for theoretical pre- 
dictions, which has been progressing steadily in recent 
years, but it also demands improvement on the observa- 
tional side with respect to defining useful observational 
samples for comparison against the best predictions that 
theory has to offer. To that end, we have undertaken a 
systematic search for giant arcs around optically selected 
clusters in two large imaging sur veys: the Second Re d Se- 
quence Cluster Survey fRCS-2 : IGilbank et al.l[20Tll) and 
the Sloan Digital Sky Survey ()York et al.ll2000[) . An ex- 
haustive visual inspection of these two surveys has pro- 
duced a combined sample of hundreds of uniformly se- 
lected giant arcs (M. D. Gladders et al. 2011, in prepa- 
ration; M. B. Bayliss et al. 2011, in preparation). Com- 
plete definition of these samples (completeness, purity, 
effective area/volume probed) will be forthcoming in fu- 
ture work, but even prior to the completion of the sample 
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Fig. 1. — The probability distribution of photometric redshifts for 
the population of foreground galaxy clusters that are responsible for 
producing the giant arcs analyzed in this paper. 



definition there is a tremendous amount to be learned by 
characterizing the properties of the giant arcs themselves. 

In this paper we use photometric color criteria to sort 
a uniformly-selected sample of 105 arcs into four red- 
shifts bins at z > 1, and evaluate the implications of 
the resulting arc redshift distribution. The redshifts and 
linear arc brightnesses of the sample are two observables 
that can be measured directly from the data, and used to 
inform future efforts to generate simulation-based predic- 
tions of giant arc abundances and match them to our new 
large giant arc samples. As discussed above, we can also 
provide direct insight regarding a possible resolution of 
the "giant arc problem" by definitively establishing that 
most giant arcs are galaxies that reside at high redshift 
(i.e. z > 2). 

This paper is organized as follows: in § [2] we summa- 
rize the photometric and spectroscopic data analyzed in 
this paper. §[3]describes our analysis methods, including 
aperture photometry of giant arcs and the color-based 
sorting of giant arcs into broad redshift bins. In § [4] we 
summarize the state of the current literature regarding 
investigations into possible contributing factors to the 
"giant arc problem" , and discuss the implications of our 
results. 

All magnitudes presented in this paper are calibrated 
relative to the AB system, via the SDSS. 

2. OBSERVATIONS 

2.1. Giant Arcs Targeted for u— band Observations 

The objects analyzed in this paper are drawn from a 
pair of comprehensive surveys for giant arcs in the RCS- 
2 and SDSS survey imaging data. As a part of a larger 
collaboration, we have systematically searched the RCS- 
2 and SDSS DR7 (Abazajia n et al.,.2009,) for incidences 
of strong lensing by galaxy clusters. The search is per- 
formed by identifying lines-of-sight in the optical survey 
data which are likely to c ontain galaxy clusters usin g 
the red-sequence algorithm (IGladders k ^[2000l[2005l) . 
The red-sequence cluster finding identifies over-densities 
of red galaxies on the sky, and provides photometric red- 
shift estimates for the galaxy clusters from the redshifted 
4000 A break. The galaxy cluster samples from the RCS- 
2 and SDSS DR7 span redshift ranges of0.2<z<1.2 
and 0.2 < z < 0.65, respectively. 
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In order to identify strong lenses, each optically se- 
lected galaxy cluster is independently inspected by mul- 
tiple experts and scored for evidence of strong lensing. 
Inspectors are also randomly made to re-score lines-of- 
sight that they have already scored, thereby providing 
data that will be used to quantify the consistency of our 
scoring, as well as our final completeness as a function 
of the averaged scores. We have also conducted exhaus- 
tive follow-up observations to quantify the purity of the 
resulting sample as a function of average score. A com- 
plete description of the SGAS and RCSGA samples are 
forthcoming (M. D. Gladders et al., in preparation, and 
M. B. Bayliss et al., in preparation, respectively). The gi- 
ant arcs used in this paper are located around 97 unique 
galaxy clusters that span the full range in photometric 
redshifts of the parent catalog of optically selected clus- 
ters described just above, with a median lens redshift of 
zi = 0.58. In Figure[T]we show the the probability distri- 
bution of foreground cluster lens redshifts for the giant 
arc sample. 

The giant arc sample presented in this paper is an in- 
complete subset of the full SDSS and RCS-2 giant arc 
samples. In the paper we analyze photometric data 
for 105 systems which constitute a large fraction of the 
final combined SGAS -I- RCSGA samples. The crite- 
ria for inclusion in our analysis here are: 1) the gi- 
ant arc was identified in a systematic search of 17— band 
imaging with depth matching the limits of the RCS- 
2 g— band survey data, and 2) that we possess pho- 
tometric u— band imaging of that arc. The u— band 
data is essential for classifying giant arcs based on well- 
established pho tmetric dropout crite ria for distant blue 
galax ies (e.g., iSteidel et al.l Il996allbl : iLowenthal et al.l 
Il997f ). All u— band observations presented in this paper 
were obtained at the 4.1m Southern Astrophysics Re- 
search (SOAR) Telescope located on Cerro Pachon in the 
Chilean Andes during the 2008B and 2009A semesters. 
In practice this means that 105 arcs discussed in the pa- 
per represent those arcs which had the appropriate a, 5 to 
be observable during the SOAR observing runs, and were 
discovered prior to the observing runs. The southerly lo- 
cation of SOAR restricts us to arcs with 5 <20 degrees, 
and because our systematic giant arc search is a pro- 
cess that is still on-going, there are many arcs which are 
now in the sample and have appropriate S to be observ- 
able from SOAR but were not known at the time of the 
u— band imaging runs. 

As discussed in the introduction, it is important to 
measure the properties for a large, unbiased sample of 
giant arcs. To that end we note that targets that were 
observed with SOAR and analyzed here were selected 
independent of their respective arc radii. Rare- The arc 
radius is measured as the mean angular distance on the 
sky between a giant arc and the center of the foreground 
cluster lens, where the center is approximated as the cen- 
ter of the BCG. The sample analyzed here includes arcs 
spanning a large range in arc radius - 3" < Rare ^ 55" 
- and the arcs were identified in imaging data that is 
approximately 1.5 — 2 magnitudes (point source) deeper 
than those used in IBavliss et al.l (|2011aD . which was re- 
stricted to arcs that were visually identified in the SDSS 
DR7 imaging data. 

Some of the giant arcs analyzed in this paper appear 
in the literature and have published spectroscopic red- 



shifts^ RCSGAJ032727- 132609 appears in lWuvts et all 
(|2010f) and iRigbv et al.l ()201lD : it is the most spectac- 
ular giant arc discovered to date in RCS-2 (and one 
of the most spectacular systems in the observable uni- 
verse). RCSGA J152745-I-G65219 is previously pub- 
lished as SGAS J15 2745+065219 (IKoester et all [20T0I: 
IBavliss et al.ll2011bh . and is located in a region where 
the RCS-2 and SDSS footprints overlap. SGAS 211119- 
011429 was fir s t iden tified as a probably giant arc in 
IHennawi et all (|200 8^. and both SGAS 095739-1-050928 
and SGAS 2 11119- 011429 have spectroscopic redshifts 
published in (Bayliss et al. 2011b). Lastly, the so-called 
"Cosmic Eye" ((Small et al.i,2007i) appears in our arc sam- 
ple, as it is located near on the sky to a massive fore- 
ground galaxy cluster, though the arc itself is formed 
around an elliptical galaxy located behind the galaxy 
cluster. We include this arc in our sample because 
it appears in our visual inspection of cluster lines-of- 
sight, and the fo reground cluster co ntributes significantly 
to the lensing (jSmail et al.l |2007[) . Several other arcs 
analyzed in this paper have been published as candi- 
date but unconfirmed giant arcs or cluster lenses, in- 
cluding SGAS 084647-1-044608, SGAS 085429-f 100819, 
SGAS 111504-H 64533, and RCSGA 004827-1-031114 
(|Wen et al.ll2"0ll . 

2.2. Follow-up u— band Imaging 

All u— band observations were conducted as part of 
NOAO programs 2008B-0400 and 2009A-0414 (PI: M. 
Bayliss) on the 4.1m SOAR Telescope with the SOAR 
Optical Imager (SOI). SOI is a small mosaic of two 
4096 X 2048 pixel CCDs filhng a field of view of 5.25 x 5.25 
arcmin, with an unbinned pixel scale of 0.0727" pixel" ^. 
Observations of all but one of our targets were taken with 
the detector binned 2x2, where the remaining target was 
observed with the detector binned 4x4. Individual ex- 
posures varied between 120s and 900s and were dithered 
to cover the chip gap in the SOI mini-mosaic. Typi- 
cal total integration times for individual targets range 
from 240s to 3780s, tuned to provide a minimum imaging 
depth complementary to the brightness of each visually 
selected giant arc in the RCS-2 g— band data, with a few 
syste ms of particular in terest being imaged more deeply 
(e.g.. lWuvts"eraIll20TOh . 

Basic image reductions were performed using a combi- 
nation of custom IDL code and the MSCRED IRAlfl 
package. The custom IDL code was employed to re- 
move time- variable detector-based noise structures that 
appeared in images, as well as for alignment of individ- 
ual exposures. Flat-fielding, bias correction, and stack- 
ing were done using standard MSCRED routines. The 
resulting stacked tt— ba nd images were calibrated from 
stars in the SDSS DR7 (jAbazaiian et al.l l2b09l that were 
observed over a range of airmasses: 1.1 < Z < 1.7. For 
each standard star observation, we solve the following 
equation: 

mu = mi,u -\- Co -\- ci X {Z - 1) + C2 X (m„_,„ J 

^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 
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TABLE 1 

Redshift Bins Definitions in Color-color 
Space 



5300 5400 5500 5600 5700 5800 5900 
K.. (A) 

Fig. 2. — Optical spectra with strong Lyman-a emission from two 
arcs with optical colors measured in this worlc that identify them as z > 
2.7 galaxies. The spectroscopic redshifts corroborate the photometric 
redshift assignment. In each panel the solid histogram represents the 
measured flux per spectral pixel, and the dotted histogram is the error 
array. Top : Optical spectrum of RCSGA 030945-143714 taken with 
the Goodman Spectrograph on the SOAR 4.1m telescope. Bottom : 
Optical spectrum of SGAS 111504-1-164528 taken with DIS on the 3.5m 
telescope at APO. 

where to„ and nig are the true apparent magnitudes of 
the standard in the u— and g— bands, respectively, m/^„ 
is the measured instrumental magnitude of the standard 
in a given u— band observation, and Z is the airmass of 
a given observation. By solving the system of equations 
resulting from all of our standard observations we deter- 
mine best-fit values for the photometric zero point cq, 
the atmospheric extinction term ci , and the photometric 
color term C2 for our u— band observations. The aver- 
age best-fit values for these terms in 2008B and 2009A 
observations is ci — 0.51 and C2 — 0.02. The limiting 
source of error in the m— band photometric zero point 
determination comes from the fundamental uncertainty 
in the SDSS w— band calibration; w— band data is noto- 
riously difficult to calibrate absolutely, a fact that has 
been well-documented in SDSS photometry 0. From the 

^ see |http://www.sdss.org/dr7/start/aboutdr7.html| 
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intrinsic scatter in our standard star observations, and 
the published photome tric characterization for the DR7 
(jAbazaiian et al.|[2009D . we adopt an upper limit on the 
uncertainty in our it— band photometry of ±0.05 mags. 
We also note that there is an established red light leak in 
the SDSS u filter that causes slight biases in photometry 
for red objects, and we therefore restrict stars used in 
our zero point calibrations to types K and bluer. 

2.3. Survey Imaging of RCS Giant Arcs (RCSGA) 
Survey Objects 

The RCS-2 survey consists of grz imaging over approx- 
imately 700 deg^ taken in queue mode with MegaCam at 
the 3.6m Canada-France-Hawaii Telescope (CFHT) be- 
tween 2005 and 2008. Individual RCS-2 pointings have 
single exposures of 240s, 480s, and 360s in each of the grz 
filters, respectively. Because these data are single point- 
ings in each filter, the images contain cosmic rays and 
chip defects which we remove manually when they occur 
near one of our giant arcs. The RCS-2 survey data is 
fully calibrated to the SDSS in all bands, which includes 
a color term describing the difference between the CFHT 
MegaCam g— band and the original SDSS 5- band. The 
RCS-2 imaging is remarkably uniform in depth and see- 
ing, making it an ideal dataset for homogeneously select- 
ing giant arcs. Fo r further detail s on th e RCS-2 data we 
refer the reader to lGilbank et al.l (|2011f) , which describes 
the RCS-2 survey data in great detail. 

2.4. Optical Imaging of Sloan Giant Arcs Survey 
(SGAS) Objects 

The primary (7— band imaging of SGAS systems was 
obtained at 2.5 — 4m class telescopes over several years 
(for examples and details of these observations see 
IHennawi et al.l [20081 ) . and it is these data from which 
the SGAS arcs analyzed here are systematically identi- 
fied. However, there exists for some of the SGAS giant 
arcs analyzed in this paper a variety of deeper imag- 
ing which we use for photometric measurements where 
it is available. A subset of the SGAS objects in this pa- 
per were also observed in the g— and r— bands with the 
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Fig. 3. — Similar to Figure 1 in lAdelberger et al.l I I2004I) . six panels sliow a synthetic galaxy spectral energy distribution (SED) at a series of 
different redshifts, with the ugrz filter transmission curves over-plotted as colored dotted lines. The ratio of galaxy fluxes measured in each of the 
different filters can reliably identify galaxies within different redshift intervals based on the location of strong spectral features, such as the Lyman- 
limit at 912A (rest-frame), and the Balmer/4000A break. These strong spectral features appear generically in star-forming galaxies and can be 
used to assign approximate redshifts for star-forming galaxies without attempting to constrain their specific stellar populations and star-formation 
histories. 



8in Gemini South Te lescope with the GMOS instrument 
(Bavhss et al.ll2011bl) and/or the 8.2m Subaru Telescope 
(Pguri ct al. 2009, 201 1]) as a part of a larger program 
designed to collect extensive follow-up observations for 
a large subset of the full SGAS sample. We use the 
best available data for all photometric measurements pre- 
sented in the paper in order to achieve the highest qual- 
ity measurements possible. Four of the SGAS objects 
analyzed here have no deep r— band imaging available 
beyond the publicly available SDSS survey imaging, and 
so we use these data where necessary. AH of the arcs, 
though confirmed in deeper band imaging, are suffi- 
ciently bright as to be well-detected in the SDSS r— band 
imaging data. 

2.5. New Spectroscopy of Individual Giant Arcs 

In this paper we also present new spectroscopic obser- 
vations of two arcs in our photometric sample. RCSGA 
030945-143714 and RCSGA 030945-14 3717 were ob- 
served with the G oodman Spectrograph ([Clemens et al.l 
12004 iCrain et al.1 [2004) on November 1, 2008 at the 
end of the first SOAR m— band observing run mentioned 
above (NO AO 2008B-0400). The Goodman Spectro- 
graph is an imaging spectrograph with multi-object capa- 
bilities designed to use Volume Phase Holographic (VPH) 
transmission gratings and optimized for throughput in 
the wavelength range, 3200 - 8OOOA, especially in the 
blue. We selected the target for spectroscopy from the 
full M— band imaging target list based on the presence 
of two bright arcs in the field (RCSGA 030945-143714 
and RCSGA 030945-143717) that could be observed si- 
multaneously by positioning the slit approximately 40 



degrees East of North so as to place both arcs within the 
slit. The spectrograph was configured with a 1.03" wide 
longslit mask, the KOSI_600 grating, and camera & grat- 
ing angles of 20 & 10 degrees, respectively. The detector 
was binned 2x2, resulting in a spatial scale along the 
slit of 0.3" pix^^ and a mean dispersion of 1.3lApix~^ 
and a spectral full width at half max (FWHM) of 4.2A. 
This mode provides a central wavelength of 5696A and 
wavelength coverage over the range, AA = 42OO-7OOOA. 

Science observations of these targets consist of 4 x 1800s 
exposures, with quartz lamp flatfields and HeAr lamp 
calibration exposures bracketing individual science expo- 
sures. The data were reduced, extracted, and stacked us- 
ing custom IDL scripts that incorporate procedures from 
the XIDL software package. At the position of RC- 
SGA 030945-143714 along the slit there is a single mildly 
asymmetric emission line feature is evident in all four in- 
dividual exposures at 4523A, which we determine to be 
Lyman-a A1216A at a redshift, Zq^c — 2.721. Other 
plausible interpretations for the emission feature (i.e. 
nebular emission lines at lower redshift) are ruled out by 
the absence of additional lines in the wavelength range 
redward of the emission feature (AA = 4600 - 7000A). 
If the emission line was, in fact, 0[II] A3727A or H- 
(3 A434lA then we would expect to also see one or more 
of H-/3 A4862A, 0[III] A4960, 5007A or H-a A6563A) in 
the wavelength range covered by our observations. We 
can also note that in our final photometric analysis RC- 
SGA 030945-143714 falls into the appropriate region of 

* http:/ /wwvir. ucolick.org/~xavicr/IDL/index. html 
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Fig. 4. — Giant arcs with known spectroscopic redshift are plotted in their corresponding color-color space. Each panel plots arcs which match 
a different redshift bin in color-color space. The bins, beginning in the upper left and proceeding clockwise are, 2.7 < z < 3.5, 1.4 < z < 2.7, 
1.4 < z < 2.1, and 1.0 < z < 1.5. There is some small overlap between some of the neighboring bins, and some arcs are therefore plotted in multiple 
panels. Arcs with known redshifts agree very well with the corresponding regions in color-color space, with 1/19 lying approximately 1.S(T outside 
of the expected region and several others lying Icr or so outside of the appropriate regions. One arc with a spectroscopic redshift, z^p^c < 1 is 
plotted as an open circle in the lower left panel, and as expected it falls outside of the color-color region corresponding to 1.0 < z < 1.5. 



the u — g vs g — r space for a galaxy with redshift z: 
2.7 < z < 3.5. At the location on the slit correspond- 
ing to RCSGA 030945-143717 there is some weak contin- 
uum signal spanning the highest-throughput wavelength 
range for the setup, but nothing sufficient to measure a 
redshift. 

A second giant arc from our sample, SGAS 
111504-1-164528, was also observed with the Dual Imag- 
ing Spectrograph (DIS) on the 3.5m Astrophysical Re- 
search Consortium (ARC) Telescope at Apache Point 
Observatory in New Mexico on the night of May 23, 
2011. These observations were conducted with DIS in 
1.5" longslit mode using the B400/R300 gratings in the 



"low-res" configuration. Science exposures consisted of 
2x2400s, with quartz lamp flatfields and HeNeAr lamp 

calibrations bracketing the science exposures. Blue-side 
spectra have a mean dispersion of SX = 1.83Apix~^ 
and a spectral FWHM of 5.3A, and red-side spectra 
have mean dispersion, S\ = 2.3lApix~^ and a spectral 
FWHM of 6.4A. The blue and red channels combine to 
provide full optical wavelength coverage over the range, 
AA = 3800 — 9800A. These data were reduced and ana- 
lyzed in the same fashion as the Goodman spectroscopy 
described above using custom IDL scripts incorporating 
procedures from the XIDL package. 
DIS uses a dichroic splitting optic centered at 
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TABLE 2 

Giant Arcs With Published Redshifts 



Giant Arc 



Reference 



RCSGA 030945-143714 2.719 

RCSGA 032727-132609 1.704 

RCSGA 213513-010143 3.07 

SGAS 095739-f050929 1.820 

SGAS 111504-f 164528 3.463 

SGAS 152745-f065219 2.760 

SGAS 211119-011432 2.858 



this work 
Rigbv et al. (2011) 
"Cosmic Eye" (Smail ct al. 2007) 
BavUss et al. (2011b) : Diehl etsT^OOg) 
tTu^wor^^^^^^ 
Koester et al. (2010) : Bavliss et al. (2011b) 
Bavliss et al. (2011b ~ 



5500A , and we observe a single strong emission line 
in both the blue- and red-side spectra corresponding to 
the slit position of SGAS 111504-1-164528, centered at 
5425A. Similar to the case of RCSGA 030945-143714 de- 
scribed above, we conclude that this feature is Lyman- 
a A1216A at a redshift, Zarc = 3.463, based on 
the lack of other visible emission features over the ob- 
served wavelength range AA = 3800 - 9800A. Similar 
to the case for RCSGA 030945-143714 above, the SGAS 
111504+164528 has u — g vs g — r colors that also iden- 
tify it as a galaxy in the redshift range, 2.7 < z < 3.5. 
Extracted spectra for RCSGA 030945-143714 and SGAS 
111504+164528 are shown in Figure [2 

3. MEASUREMENTS AND METHODOLOGY 

3.1. Optical Selection of Star Forming Galaxies 

Historically, Lyman Break Galaxies (LBGs) are se- 
lected by examining wide-band photometric data and 
identifying the redshifted 'Lyman limit' continuum 
break. This str ong spectral break a pp ears at 912A in 
the r est frame (jSteidel et al.l Il996al lbl: iLowenthal et al.l 
[1993), and moves redward in the rest-frame with increas- 
ing redshift - approaching 1216A - due to Lyman-a for- 
est (jSteidel fc Sargent! [l987l : iRauch 1998) absorption by 
intergalactic neutral hydrogen. Surveys for LBGs are ef- 
ficient for collecting statistical samples of high-redshift 
galaxies because it is difficult for galaxies at lower red- 
shifts to mimic the s udden and extreme spectral break 
of the Lyman Limit. iSteidel et al.l ([2003) published the 
comprehensive results of a systematic search for z > 2.7 
galaxies, and this work was extended to probabilisti- 
cally sort star-forming galaxies into bins at lower red- 
shift intervals by identifying the color-evolution that oc- 
curs as additional weaker spectral features redward of 
the Lyman Limit redsh i ft through a set of opt ical fil- 
ters ([Steidilitin llOOl |AdeIberge^^ The 
ratio of galaxy fluxes - or colors - measured in differ- 
ent broadband near-UV, optical, and near-infrared filters 
can reliably identify galaxies within different redshift in- 
tervals based on the location of strong spectral features, 
such as the Lyman-limit at 912A (rest-frame), and the 
Balmer/4000A break (see Figure [3]). These strong spec- 
tral features appear generically in star-forming galaxies 
and can be used to assign approximate redshifts for star- 
forming galaxies without attempting to constrain their 
specific stellar populations and star-formation histories. 
We use this color selection technique and use it to iden- 
tify the fraction of a complete sample of 105 giant arcs 
which fall into four broad redshift bins at z > 1 (and by 
exclusion, giant arcs at z < 1 are also identified). 



In order to defin e photometric dropo u ts we follow 
the rnethod ologv of ISteidel et all (I2003D . ISteidel et afl 
(|200l . and lAdelberger et al.l (|2004[ )~to define regions 
in color-color space that correspond to distinct redshift 
bins. The available photometric data i s in a set of filter s 
that are similar to those used in Steidel et al. (2003), 
St eidel et al. (2004), and Adelberger et al. (2004), with 
the only significant differences being an r— band filter 
with somewhat different ch aracteristics. We u s e syn - 
thetic galaxy spectra from iBruzual fc Chariot (|2003( ) 
with identical stellar age and star formation history prop- 
erties as are used by Steidel et al. (2003), Steidel et al.l 
(|2004f ). and lAdelberger et all (120041 ) and step these 
model spectra through redshift steps of Sz = 0.01 
from z = 0.5 to z = 4.0 and measure the resulting 
synthetic colors. These simulated galaxy colors as a 
function of redshift are e s sentia ll y identical t o thos e 
appearing in 'Steidel et al." C 2003D . Isteidel et all (|2004l ). 
and [Adelberger et al. (2004=) with slight differences of 
^ 0.1 — 0.15 mags in colors which include the r— band fil- 
ter. The final color cut criteria that we apply to our pho- 
tometry defines four primary redshift bins, all of which 
are summarized in Table [TJ 

Arcs that do not satisfy any of the criteria in Table [T] 
are likely to have redshifts that are either greater than 
3.5 or less than 1.0. In principle we might be able to iden- 
tify very high redshift arcs in our sample (e.g. z > 3.5) 
as g— band dropouts using logic that is ana logous to the 
origin al m— band dropout definitions from ISteidel et al.l 
((200l . However, we are limited by the lack of deep imag- 
ing redward of the r— band. In any event, arcs that are 
very red in g — r are rare {Ng-r>i.2 — 6), and all of these 
objects either have r — z > 0.9, or are detected in z at less 
than 4cr. This leads us to conclude that arcs which are 
red in g — r are likely to be intrinsically red galaxies at z 
< 1. That stated, it is still possible that our arc sample 
contains one or two very high redshift g— band dropout 
galaxies, and we account for this possibility in the final 
uncertainties on our constraints for the giant arc redshift 
distribution. 

3.2. Aperture Photometry of Giant Arcs 

We measure four- and three-band photometry for a full 
sample of 105 giant arcs from the RCSGA and SGAS 
samples, respectively. All objects have measurements in 
ugr, and objects located in the RCS-2 fields have z— band 
measurements as well. Photometric measurements are 
conducted using a custom IDL pipeline that allows us to 
construct photometric apertures which follow the ridge- 
lines of giant arcs and match the highly variable struc- 
tures of these objects on the sky, following the proce- 
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dure used by iWuvts et alj ()201Cll ) . We draw apertures 
by manually defining a ridge-line along each giant arc, 
and convolve the ridge-line with the point spread func- 
tion (PSF) of each image to create a series of apertures 
of increasing radius, following isophotes of the convolu- 
tion. Arcs are typically only resolved along a single axis, 
because the lensing magnification acts almost exclusively 
along the tangential direction (relative to the center of 
the potential of the foreground lens) . We therefore mea- 
sure photometry in each image out to an aperture radius 
equal to 2.0xFWHM, where the FWHM describes the 
PSF measured from reference stars in each image. These 
apertures allow us to measure an equivalent region on 
the sky for images with different PSFs. All magnitudes 
are aperture corrected to a radius of 6" using the curve 
of growth for reference stars in each image. 

In cases where the flux for a given arc within its 
2.0xFWHM defined aperture is less than 2aapp, where 
Capp is the noise measured within the aperture, we mea- 
sure a limiting magnitude for that object in that fil- 
ter. Limiting magnitudes are measured as the magnitude 
which corresponds to 2aapp, and this limiting magnitude 
is aperture corrected to the same 6" radius as above. 
Quoted limiting magnitudes are 95% confidence upper 
limits on the flux from the corresponding arc within its 
apert ure. All photometry is corrected for galactic extinc- 
tion (jSchlegel et al.lll998[ ). 

The ridge-line apertures used for this analysis are not 
necessarily intended to enclose all of the observable flux 
from each lensed background source. Rather, we draw 
apertures along the largest contiguous image or arc of 
each lensed galaxy, and make no attempt to add flux 
from multiple images/arcs, or faint arc tails that appear 
to originate or extend away from the same source. Our 
measurements are conducted with robust color estimates 
as our primary goal, and our ground-based imaging typ- 
ically does not have either the depth or image quality 
to unambiguously identify complete arc families in the 
same way as can be achieved with adaptive optics or 
HST quality imaging. 

We also note that variations in color that can be 
observed in different lensed images of a single back- 
ground source may pr ovide a distorted (e.g. cB58; 
iWilliams fc Lewii 119961) representation of the intrinsic 
source. That stated, color gradients for star forming 
galaxies te nd to be small - on the order of ~ 0.1 mags 
(jSuh et aT. 2010). Color distortions due to magnification 
gradients across the surface of the lensed sources should 
therefore be well within the color variations accounted 
for in the definition of the regions in color-color space, 
which are intended to apply to star forming galaxies with 
a variety of colors and there is no reason to expect catas- 
trophic outliers (i.e. giant arcs with apparent colors that 
differ dramatically from their true intrinsic colors) to be 
a significant problem for our analysis. 

Approximately 1 / 3rd of the arcs analyzed here are lo- 
cated nearby on the sky to other, presumably foreground, 
galaxies. For these objects there is light from a projected 
source, typically the brightest cluster galaxy (BCG) or a 
cluster member galaxy in the foreground lens, that over- 
laps with the 2.0 X FWHM aperture that we use to mea- 
sure giant arc arc magnitudes. We have used the pub- 
licly available GALFIT package ()Peng et al.ll2010D to fit 
Sersic profiles to the nearby galaxies and subtract the 
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Fig. 5. — Photometric colors, u — g vs g — r, for all 105 giant arcs that 
have It — band data from SOAR. Three regions in the color-color space 
are defined by the dashed, dotted, and dot-dashed lines corresponding 
to rcdshift ranges, 2.7 < z < 3.5, 1.9 < z < 2.7, and 1.4 < z < 2.1, 
respectively. Arcs that fall into each respective bin arc plotted as blue, 
cyan, and green points, respectively, while arcs not falling into any 
of the aforementioned rcdshift bins arc plotted in black; by exclusion 
these objects have redshifts < 1.4. 



resulting model fluxes. Galaxy light profiles are fit in 
a filter where the galaxy signal-to- noise (S/N) is high, 
and that profile is then saved and scaled in brightness 
to match the galaxy light in the other imaging bands. 
A single Sersic profile is almost always sufficient to the 
task for this work because we are generally concerned 
with subtracting off the extended stellar light profile of 
BCGs, and we are not concerned with accurately mod- 
eling the detailed structure in the BCG core which has 
no impact on the arc photometry. Occasionally there are 
very faint intervening galaxies that lack sufficient S /N to 
achieve a robust GALFIT model, and we remove these 
galaxies by manually masking the galaxy pixels. 

3.3. Verification of Color-Based Redshift Assignment 

We can directly test the validity of ap plying 
the iSteidel et all (l200l . iSteidel et all (pOOl . and 
lAdelberger et al.l (|2004[ ) color criteria for rough redshift 
assignment to our sample of giant arcs by comparing 
spectroscopic redshifts for individual giant arcs against 
their color-color space redshift designation. All photo- 
metric measurements were made blind to any prior red- 
shift knowledge for our giant arc sample, and after the 
fact we identify 19 giant arcs in our sample that have 
known, high-confidence spectroscopic redshifts. Five of 
these giant arc redshifts are in the literature, while the 
remaining 14 are the result of new observations. Two of 
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Fig. 6. — Colors, g — r vs r — 2, for 96 arcs discovered in the RCS-2 
survey area. The dashed line indicates the division in the color-color 
space that identifies typical star forming galaxies at 1.0 < z < 1.5 
(these galaxies are located to the right/below the dashed line). Arcs 
satisfying this criteria are plotted in purple, while other arcs arc plotted 
in black. Several arcs, including many of those with marginal or poor 
2— band photometry fall outside of the color-color region shown in the 
plot. 



these new redshifts are presented here (see Section 2.5), 
and the remaining 13 will appear in forthcoming publica- 
tions: one in M. D. Gladders et al. 2011, in preparation 
and 12 in M. Carrasco et al. 2012, in preparation). The 
authors of these papers have generously granted us ac- 
cess to the unpublished spectroscopic redshifts for the 
purpose of verifying the color-color redshift designation 
method used in this paper. The seven giant arcs with 
public redshift information (five published plus two new 
arcs presented here) are summarized in Table [2j 

In order to verify the success of the color-color based 
methodology when applied to our giant arc sample, 
we populate the relevant color-color spaces with mea- 
sured colors for arcs of known redshift in our sample 
and compare their location to the regions define d base d 
on the meth o ds fir st pu blished in iSteidel et al.l (I2003D . 
ISteidel et all (|200l . and lAdelberger et al.l ()200l . The 
resulting u — g vs g — r and g — r \s, r — z plots are 
shown in Figure 21 and provide excellent validation of 
the use of cuts in color-color space as a tool for estimat- 
ing the redshifts of our giant arc sample. We can perform 
an additional sanity check by comparing the redshift bin 
assignments for the giant arcs against photometric red- 
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shifts for the foreground cluster lenses that are measured 
as a part of the red-sequence cluster finding algorithm. 
Indeed, we find that the giant arcs which are sorted into 
the z < 1 redshift bin are not formed by clusters that 
have Zphot ^ 0.8, which is encouraging given that the ge- 
ometry of such systems makes the occurrence of strong 
lensing difficult at best (and non-physical at worst). A 
broad examination of the relationship between ziens/^arc 
pairs yields no significant correlation. 

3.4. Redshift Constraints On 105 Giant Arcs 

Giant arc u — g vs g ~ r and g — r ys r — z colors 
are shown in Figure [5] and Figure [51 respectively. Indi- 
vidual giant arc color measurements are color-coded in 
the figures according to the redshift bins that they fall 
into. The g — r vs r — z color space in Figure [S] is popu- 
lated using only RCSGA survey objects, because we lack 
z— band imaging of sufficient depth for our SGAS arcs. 
RCSGA arcs, however, comprise > 90% of our full giant 
arc sample, and we can also note that seven of the nine 
SGAS arcs fall into redshift bins based on their u — g 
vs g — r colors, so that the lack of z— band data for the 
SGAS sample only prevents us from assigning two sys- 
tems to one of the z < 1.0 or 1.0 < z < 1.5 bins, which 
has a negligible affect on our ability to characterize the 
redshift distribution of the sample as a whole. 

The resulting constraints on the fraction of arcs falling 
within each redshift bin are shown in Table |31 We es- 
timate the errors on the fractions of giant arcs meeting 
various redshift criteria by generating IC* Monte Carlo 
realizations of our data, where each realization is a sim- 
ulated dataset that is generated by assuming that each 
photometric measurement is gaussian distributed about 
the measured value, with standard deviations equal to 
the 1(7 errors on the measurements. Some photomet- 
ric measurements have asymmetric errors, and for these 
cases the simulated data accounts for the asymmetry in 
the errors. We make no attempt to recover a total prob- 
ability distribution, dpa,T:c/dz, because doing so requires 
detailed knowledge of the probability distribution of red- 
shifts within each bin in color-color space, which we do 
not have. Instead we simply report the total cumulative 
fraction of giant arcs which have redshifts, 7,arc^ greater 
than various discrete cuts. The resulting cumulative red- 
shift fractions are 7lt\ % at z > 1.0, 64+4 % at z > 1.4, 
56t^ % at z > 1.9, and 2ll^ % at z > 2.7. The fraction 
of arcs with z> 1.4 and z> 1.9 requires making some 
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assumptions about the probability distribution within 
each redshift bin, because some of the neighboring bins 
have regions of overlap in redshift. For example, in or- 
der to compute the fraction of giant arcs with z> 1.4, 
we first account for all arcs that have colors identifying 
them as galaxies at 2.7 < T.arc < 3.5, 1.9 < Zarc < 2.7, 
and 1.4 < Zarc < 2.1. We then must estimate the frac- 
tion of arcs within the redshift range, 1.0 < Zarc < 1-5, 
that we expect to have Zarc > 1-4. For these reported 
values we assume that giant arc redshifts within a bin 
are roughly gaussian distributed within that bin with a 
standard deviation that is based on the results of spec- 
troscopic follow-up for gala xies falling within these same 
color-based redshift bins in lAdelberger et all ()2004D . 

3.5. Giant Arc Linear Brightness Measurements 

In addition to the optical colors of our giant arc sample, 
we also measure the linear brightnesses for each arc in the 
g— and r— bands. The linear brightness is computed for 
a given arc by dividing the total aperture magnitude for 
each arc by the length of the aperture, in arcseconds. 
This quantity is a very useful observable for giant arcs 
with ground-based imaging data, where the arcs are typ- 
ically unresolved in the radial direction (with respect to 
the center of the foreground lensing potential) . So long 
as an arc is resolved along the tangential direction (again, 
with respect to the center of the foreground lensing po- 
tential) - a condition that is essentially always satisfied 
- then the linear brightness is an observable that can 
be extracted from the data. The linear brightness has 
the appealing property of being independent of the PSF 
of the seeing of a particular image, and it is an observ- 
able quantity that can be easily produced for simulated 
giant arcs by convolving simulated imaging data with 
a smoothing kernel that is chosen to match the image 
quality of real imaging data. 

Figure [7] shows the distribution of both integrated and 
linear magnitudes for our complete arc sample in both 
the g— and r— bands. The measured distribution of lin- 
ear magnitudes for a statistical sample of giant arcs is 
an observable which, like the arc redshift distribution, 
informs predictions of giant arc statistics by providing 
direct constraints on the photometric properties of sim- 
ulated giant arcs. Specifically, the linear magnitudes of 
arcs should depend on some combination of two factors: 
1) the typical magnifications for the lens-source systems, 
and 2) the intrinsic surface brightness function of the 
lensed background sources. Simulated arc samples must 
be able to reproduce the observed linear magnitude dis- 
tribution in order to ensure that they are providing a 
robust simulation of the physical universe. 

4. DISCUSSION AND SUMMARY 

4.1. A Census of Physical Factors Relevant to 
Predictions for Giant Arcs 

As we summarized in the Introduction, there has been 
no shortage of ideas for physical effects which, if un- 
accounted for in cosmological simulations, can explain 
varying degrees of discord between giant arc statistics 
predictions and observations. One such factor that has 
been directly addressed in the literature is the impact 
that baryons can have on the strong lensing c ross-section 
for galaxy clusters (|Puchwcin et al.. ,2005t iRozo et all 



[2008t iWambsganss et all [2001 . iRozo et all (|2008l ) find 
that the inclusion of baryonic cooling results in signifi- 
cant steepening of the central density profile relative to 
what is seen in dark matter only simulations and estimate 
that these processes can increase giant arc abundances 
by as much as a factor of a few, and iWambsganss et al.l 
(12001 find more modest increases of ~ 25%. However, 
the real-world impact is very likely to be smaller than 
measured by these authors due to the w ell-established 
"over -cooling" problem in simulations (e.g..lBalog h et aP 
|2001| ). In fact, more recent studies of strong lensing effi- 
ciencies in simulations that include gas physics find that 
the inclusion of feedback processes serves to mitigate the 
increase in strong lensing efficien cy that result from in- 
cluding baryonic cooling physics ()Mead et al.ll2010l) . 

One key property of massive halos that does have 
a large impact on their ensemble lensing efficiencies is 
the triaxiality of their matter distributions. Halos iden- 
tified in dark m atter simulations are k nown to have 
triaxial shapes (jKasun fc EvrardI [20051: lAllgood et all 
|2006[) . and the triaxiality can cause large orientation- 
dependent variations in the strong lensing cross-section 
for individual halos ( Oguri ct al. 2003; Da]al ct al. 200j; 
iHennawi et aD [20071 : iMeneghetti et al.H2010( ). This ef- 
fect should be well-accounted for in predictions of gi- 
ant arc counts usin g large volume s imulations (e.g., 
IHennawi et"al]|2007l: IMeneghetti et al.ll2010D . We must 
also point out, however, that the triaxiality of the mat- 
ter distribution in the cores of clusters in simulations is 
sensitive to the presence/absence of gas physics. Bary- 
onic cooling processes cause the shape of the density pro- 
file in cluster cores to become much less triaxial than in 
pure dark matte r simulations (jKazantzidis et all 120041 : 
IRozo et all 120081 ). which could mitigate the importance 
of triaxiality and the magnitude of an orientation bias. 

Simulations that account for the presence of baryonic 
matter by populating dark matter simulated halos with 
evolved mock galaxies may also be able to increase the 
number of expected giant arcs, but different studies claim 
varying increases in the efficiency of giant arc production. 
Dark matter only simulations with galaxies "painted on" 
after the fact cannot account for the detailed effects that 
baryons have on the shape and steepness of the matter 
density profile, but the inclusion of a large concentration 
of baryonic matter in the cores of simulated halos does 
produce an increase in the strong lensing cross-section. 
The size of the impact that painted on galaxies have on 
giant arc counts is claimed to be as sniall as < 25% 
(jFlores et al][2000l: IMeneghetti et al.ll2000l . [2001 . and as 
large as a factor of 2 ( Hilbert et al.ll200'8 ). Though the 
precise magnitude of the effect of including galaxies in 
ray-tracing simulations of giant arcs varies in the litera- 
ture, there is no evidence that it approaches the factors of 
~ 10 necessary to account for the diffe rence between cur- 
rent a rc counts and the predictions bv lBartelmann et al.l 
([1991 . 

The impact of dark matter substructure in cluster-scale 
halos has also been suggested as a potentially important 
factor (Fl ores et al.l 2000: Mcncghctti ct al. 2000, 2003), 
but work done to quantify the contribution of substruc- 
ture to lensing effici ency suggests that it is negligible 
(jHennawi et al.|[2007D . Other recent work indicates that 
the efficiency for giant arc production of individual clus- 
ters seems to be in good agreement when comparing the 
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Fig. 7. — Distribution of magnitudes for all giant arcs presented in this paper. Left : Integrated magnitudes for all ares; the solid line histogram 
shows g — band magnitudes, and the dashed line histogram shows r— band magnitudes. Right : Linear magnitudes for the entire are sample, where 
the solid and dashed line histograms again eorrespond to the g— and r— band, respectively. 



number of giant arcs produced per cluster in HST imag- 
ing of X-ray selected cluster samples against mock images 
creat ed from ray tracing mode rn simulated massive clus- 
t ers (IHoresh et alJ[2nn5l l2?)TTh. 
iHoresh et al.~(l201lD use the Millennium Simulation 



(jSpringel et al 



I2005D with galaxies included via a semi- 



analytic prescription, so that the galaxy contribution to 
the tensing efficiency of individual clusters is accounted 
for. Notably, this work does not address the question 
of the overall abundance of giant arcs - e.g. giant arc 
counts over a well-quantified fraction of the sky imaged 
to a well-defined depth - but it does isolate the question 
of whether there is a systematic discrepancy between the 
strong tensing efficiencies of individual simulated vs. real 
clusters. As summarized above, a study of the global 
increase in giant arc counts in the Millennium Simula- 
tion, including galaxies, suggests a factor of 2 in crease 
over the earl i er wor k of iBartelmann et al.l ()1998D . The 
IHoresh et al.l (|2011[ ) result therefore imply that the ap- 
parent discrepancy between the real and predicted giant 
arc abundances is primarily due to factors other than 
those having to do with the internal structure of galaxy 
clusters (e.g., contributions to the project surface mass 
density from abnormally high concentrations, substruc- 
ture, or cluster galaxies). 

With clear motivation now to focus on factors beyond 
the internal properties of simulated clusters, there are 
several potential effects on which to focus our attention. 
Two additional ideas that have been explored in the lit- 
erature are the contribution to cluster tensing efficiency 
from large scale structure (e.g. filaments and uncorre- 
lated halos) along the line of sight looking out toward 
mass i ve clusters dWambsg anss et al. 2005; Hilbcrt ct al. 
I2007t iPuchwein fc Hilbertll2609l ). and short timescale in- 
creases in the tensing cross-s ections of i i idivid ual clus- 
ters due to major mergers (iTorri et al.l l200l . Stud- 



ies using the Millennium Simulation conclude that the 
resulting increase in the strong tensing efficiency from 
these unrelated, projected density effects is in the 30% 
range, and are therefore not sufficient to significantly 
ease the tension between giant counts observed and pre- 
dicted. Observational tests looking fo r projected line-of- 
sight structure towar d strong lenses ()Faure et al.ll2009l : 
iFassnacht et al.l[2011 ^ find no evidence of an increase in 
the amount of uncorrelated structure near lenses vs. the 
field, in agreement with th e simulation-based results. 

With regard to mergers, iTorri et"an ()2004f ) find an in- 
crease of nearly an order of magnitude in an individual 
cluster's strong tensing cross-section during a merger, the 
short timescale of the boosted cross-section and the in- 
frequency of major mergers in cluster-scale halos make 
it difficult for mergers to contribute substantially to 
the integrated efficien c y of galaxy cluster strong tensing 
(|Hennawi et al.ll2007D . iFedeli et all ()2006D conclude that 
accounting for mergers in simulations can nearly double 
the strong tensing optical depth for a subset of lenses - 
those at z > 0.5. However, a large fraction of cluster 
lenses lie at z < 0.5, thus reducing the total magnitude 
of merger effects on simulated giant arc samples. 

Having accounted for an exhaustive set of physical 
factors that could increase giant arc counts by making 
simulated clusters into better strong lenses, we are left 
with one hypothesis that h as nothing at all to do with 
the cluster lens population. iHamana fc Futamasd (|1997l ) 
first argued that giant arc statistics depend sensitively 
on any evolution of the galaxy lun iinosity function at 
high redshift, and lOguri et at.! (|2003f) identify the uncer- 
tainty in the galaxy luminosity function at high redshift 
as a factor that can significantly impact predictions for 
giant arcs based on their semi-analytic models. This ef- 
fect was quantified in simulations by placing background 
sources at a series of different redshift planes and com- 
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puting the integrated strong lensing cross section for a 
simulate d cosmological volume as a function of so urce 
redshift (|Dalal et al."2004VWambsganss eLaI||200j). 

Interestingly, Wambsganss ct al. (2004 ) are able t o ap- 
proximately reproduce the Bartclnia nn et alJ (1998) pre- 
dictions by using a background galaxy population placed 
entirely at Zs — 1, but when background galaxies are 
placed at Zj = 1.5 they recover a factor of ~ 10 increase 
in the number of predicted giant ar cs, and a fa c tor o f 
~ 20 increase for sources at z= = 2. IDalal et alJ (l200l 
also identify a significant increase in the total strong 
lensing cross-section with increasing source redshift, but 
measure the increased cross-section a t higher source red- 
shift t o be a factor of ~ 3 smaller than lWambsganss et alJ 
(|2004f ). Simulations which evaluate lensing efficiency 
with all sources placed at a single s ource plane are not 
physically r ealistic, but the work by 'W ambsganss et aH 
HqM) and IDalal et all (|200l emphas izes the impor- 
tance of understanding the distribution of source red- 
shifts of background galaxies available to be strongly 
lensed by foreground galaxy clusters. 

4.2. These Results In the Context of Other Work 

There is strong evidence that the total strong lensing 
cross-section for galaxy clusters can easily be increased 
by as much as a factor of at least ^ 7 — 10 if a majority 
of the sources lensed into arcs reside at z > 2. Our re- 
sults analyzing the optical colors for a large, well-selected 
sample of 105 giant arcs provide a measurement of the 
distribution of redshifts for sources which are lensed into 
arcs, with 7ltl % at z > 1.0, 64^^ % at z > 1.4, 56^1 % 
at z > 1.9, and 21^2 % at z > 2.7. The inferred median 
redshift here is approximately Zs = 2.0, which is consis- 
tent with the median reds hift of 1.82 determine d from 
a much smaller sample in iBayliss et al.l (|2011al ). Also 
encouragingly, the fraction of gia nt arcs measured at z 
> [1.0, 1.4, 1.9] in this paper and in lBavliss et al.l (|2011aD 
are consistent to within the statistical errors of the two 
measured redshift distributions. As a caveat to the sort- 
ing of arcs into precise redshift bins, we point out that the 
assumptions that go into defining regions in color space 
that associate with specific redshift ranges are likely to 
be much more robust in some cases than others. Specif- 
ically, the Lyman-a forest and the Lyman-limit break 
moving through the u— band at z > 1.9 will strongly af- 
fect the u—g color of any galaxy regardless of the shape of 
its spectral energy distribution (SED). Color criteria at 
z < 1.9, however, rely more heavily on the robustness of 
the assumptions made about the properties of the stellar 
populations within the observed galaxies. 

With robust measurements of the giant arc redshift 
distribution for two samples of giant arcs, it is also im- 
portant to confirm that the measured distributions are 
sensible in the context of other observat i ons of high red- 
shift galaxy populations. iBavliss et al.l (|2011aD test the 
redshift distribution of 28 bright giant arcs against sim- 
ple physical models based on results from strong lensing 
simul ations and galax y catalogs from the COSMOS sur- 
vey ( Ilbert et al.ll2009l ) and find good agreement with the 
observed redshift distribution. The form of the model 
redshift distribution is, 

, dn 

"Pare _ "arc^^s 

dz, JdZsa^rcj^ ' 



where the mo del uses a normali zed strong lensing 
cross-section from lFedeli et al.l (|2010D evaluated at a sin- 
gle source plane, Zg = 2, and we solve for the giant arc 
cross-section at different source redshifts by assuming a 
universal density profile that is described by a power- 
law slope, a, in the regions where strong lensing is ob- 
served. The number density of background galaxies as 
a function of redshift can be computed from the COS- 
MOS catalogs for a given limiting magnitude and an es- 
timate of the average magnification factor. We do not 
have good constraints on the average magnification of 
these sources, and the literature is not unanimous on the 
precise value of the slope of the density profile, a, in the 
cores of clusters. However, by testing models that sam- 
ple a range of reasonable values for these two parameters 
(—1.7 < a < —1.3, 26 < giim < 24) we can produce 
synthetic redshift distributi ons that are in good agree- 
ment with the data (Bavliss et al.|[2011aD . Though the 
broad bins that we define in this paper limit our ability to 
pe rform the same statis tical comparisons that are used 
m IBavliss et all (|2011aD . we can compare the fractions 
of arcs above various thresh-holds in the models and in 
our color-based redshift constraints and find good agree- 
ment. Because the sample analyzed in this paper was 
identified in imaging that is approximately 2 m agnitudes 
deeper than that used in lBavliss et al.l ()2011aD . the fact 
that both redshift distributions can be well-described by 
the same simple physical models has an interesting impli- 
cation; the primary difference between the brightest arcs 
selected from shallower imaging vs somewhat fainter arcs 
found in deeper imaging is that the brighter arcs are typ- 
ically the cases with the most extreme magnifications. 

4.3. Implications for Giant Arc Statistics 

With the giant arc redshift distribution now con- 
strained in this paper, it is possible to calculate the 
approximate increase in the total cross-section for arc 
production compared to published predictions. We 
estimate the relative cross-section for different back- 
ground source redshift distributions using Figure 1 in 
IWambsganss et a l. (2004i), which provides a measure- 
ment of the strong lensing cross-section as a func- 
tion of source redshift. The total cross-section for our 
measured redshift distribution is simp l y the integrated 
cross-sections from IWambsganss et al.l (|2004D Figure 1, 
weighted by the fraction of giant arcs determined to be 
within different redshift bins. The uncertainty in this de- 
termination is propagated from the color-based redshift 
identification, estimated by varying the weights in accor- 
dance with our measured uncertainties in the fraction of 
giant arcs located within each redshift bin (see Table [S]) . 

The resulting total increase in the efficie ncy for gi- 
ant ar c production, normalized to the Bartclm ann et ahl 
(|1998f ) result with lensed sources at Zs = 1, is 
(tsl/ctsl, 2=1 = 10.3^Qg. The uncertainty reported in 
this measurement does not incorporate an attempt to 
account for published evidence that the magnitude of 
the effect i s smaller than Wambsganss et al. (2 004) claim 
(|Dalal et al . 2004; Li ct al. 2005rFcd eh et al.ir20G6l ). but 
it is clear that the high-redshift nature of the lensed 
sources is a significant factor to be accounted for in 
predictions of giant arc abundances. We can compare 
our new estimate of <Jsl/o'sl,z=i to the same value es- 
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Fig. 8. — The probability distribution for tlic estimated increase 
in the total cross-section for the production of gravitationally Icnsed 
arcs with length-to-width ratios > 5, relative to the predictions 
of Bartclmann ct al.l <|1998I ). Estimates are made from the redshift 
distribution derived in this paper using color-based redshift assign- 
ments (solid line), as well as the redshift distribution reported in 
IBavliss et al.l ]2011al) (dashed line). The two estimates of the increase 
in giant arc cross-section have nearly identical expectation values of 
(y SL / <y SL ,z=l — 10.3 — indicated by the vertical red line — but the 
much larger sample used in this work vastly improve the precision of 
the estimate. 



timated using the spectroscopic redshift distribution of 
IBavliss et ahl (pfiTTal ). (JsL/(rsL,z=i = lltg (Figure |8|. 
The improved constraints on the arc redshift distribu- 
tion using a sample of 105 - along with the fact that 
four of the 28 arcs analyzed in IBavliss et al.l (|2011aD had 
only very broad redshift constraints spanning a range in 
Az = 1.8 - provides a much tighter estimate of the in- 
crease that we should expect in giant arc counts when 
using an empirically calibrated redshi ft distribution, rel- 
ative to the simple Zg = 1 used in iBartelmann et ahl 
()1998[ ). Modern efforts to simulate giant arc counts can 
incorporate the empirically determined giant arc redshift 
and brightness measurements that we provide in this 
work, and thus produce refined predictions for giant arc 
statistics. 

We have established that incorporating a background 
source population into future efforts to predict giant 
arc counts should increase the number of predicted arcs 
by as much as a factor of ~ 10 beyond the results of 
IBartelma nn et aH (1998), and thereby diffuse the ten- 
sion between observations and prediction for arc counts. 
However, we have yet to discuss an important detail re- 
garding most published work in which giant arcs are pro- 
duced by ray-tracing simulations. Specifically, most sim- 
ulations run in the past ~ 10 — 15 years were run in 
cosmologies with a value for the normalization of the 
matter power spectrum, as, that is higher - typically 
(Tg = 0.9 — 0.95 - than the value that is preferred by 
current observational constraints of erg 0.81 ± 0.03 
(|Komatsu et al.l 120 lit iLarson et~aIll201lD . The precise 
value of ug used in simulations is of particular impor- 



tance for predictions of giant arc statistics because the 
high-end tail of the mass function varies sensitively with 

The quantitative impact of erg 0.81 on predi ctions of 
giant arc counts has been studied wit h N-body (Li et alj 
120061 ) and semi-analytic techniques ([Fedeli et al . 2008|) , 
and both confirm the intuitive expectation that lower 
values of erg result in significantly lower giant arc abun- 
dances. The difference between a cosmology with erg — 
0.9 vs (7g = 0.7 4 is found to be a factor of 6 in sim- 
ulations (ILi et a l. 2006), and semi-analytic modeling by 
IFedeli et al.l ()2008D suggests that tjg ^ 0.8 under-predicts 
giant arc counts by an order of magnitude even after 
attempting to account for factors such as mergers and 
cluster su bstructur e. 

The Fed eli et al.l (j2008f) models do not account for con- 
tributions from galaxies embedded in cluster potentials, 
and use a background source description that is based on 
somew hat out-dated observat ions using the Hubble Deep 
Field (jCasertano et al.l I2000D that are subject to large 
cosmic variance uncertainties (jHennawi et al.|[2007[ ) , but 
the fact that an order of magnitude discrepancy persists 
in this work is clear motivation for the pursuit of mod- 
ern simulations to better-characterize the expected giant 
arc abundance. The machinery is in place to produce 
these predictions in such a way that any persistent under- 
abundance of giant arcs in theory can be reasonably in- 
terpreted as evidence for real tension between the growth 
of structure in ACDM- in the form of the abundance of 
the most massive clusters - and the integrated strong 
lensing properties observed in the universe. New, large 
giant arc samples are in the process of being identified 
from wide-field imaging survey data, and will soon be 
available in the literature (M. D. Gladders et al. 2011, 
in preparation; M. B. Bayliss et al. 2012, in prepara- 
tion). These forthcoming samples will provide observa- 
tional constraints that match the fidelity of the next gen- 
eration of improved predictions, and set the stage for a 
definitive test of giant arc statistics. 

4.4. Other Applications for A Catalog of Lensed 
Galaxies At z > 2 

The current literature contains a modest sa mple of ex- 
treme ly bright strongly lensed galaxies (e . g.. lYee et al 
199 6; Small et al.ll2007[|Koester et al.|[2010HWuvts et al 



|2010|), some of which have been the target of d etailed 
multi- wavelength follow-up studies fiFin kelstei n et al.l 
I2OOI ISiana et all [20091 : IQuider et all [201(1 . Studies of 
individual galaxies at high redshift are limited by the in- 
trinsic faintness of the sources - especially where high fi- 
delity observations of typical ^L* galaxies are concerned. 
However, because faint galaxies (L<L*) are far more nu- 
merous than bright ones and therefore more likely to lie 
on the caustic of foreground lensing clusters, we should 
expect most giant arcs to be intrinsically sub-L* galaxies. 
It follows then that strongly lensed galaxies are ideal tar- 
gets for studying the properties of typical high-z galaxies 
at S/N that is unavailable in studies of galaxies in the 
field. By analyzing the optical broadband colors of 105 
arcs, we have shown that approximately half of all gi- 
ant arcs formed around galaxy clusters can be expected 
to lie at z > 2. Forthcoming large catalogs of hundreds 
of galaxies that are strongly lensed by clusters should 
therefore provide of order hundreds of strongly lensed 
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high-redshift sources that will be accessible for detailed 
individual study with the instrumentation available on 
8 — 10m class telescopes. 



4.5. Summary 

We have measured ugr magnitudes for a sample of 105 
giant arcs, and z— band magnitudes for 96/105 of the 
same sample. The giant arc sample was identified in 
a systematic search of uniform depth g— band imaging 
data, and constitutes a representative subsamplc of two 
forthcoming large samples of giant arcs that are currently 
being identified in search of well-characterized cosmolog- 
ical search volumes. We then sort all 105 arcs into five 
redshift bins based on well-established broadband color 
criteria, and thereby constraint the underlying redshift 
distribution of giant arcs identified with our selection 
method. The data indicate that Zg = 2.0 ±0.1 is the typ- 
ical redshift for these giant arcs, which agrees remark- 
ably well with the redshift distribution measured from 
spectroscopic observations of a much smaller sample (N 
= 28) giant arcs that were identified in much shallower 
imaging data. By establishing that approximately half 
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